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ABSTRACT 

Using the Keck Interferometer, we spatially resolved the orbit of the 
pre— main-sequence binary, Haro l-14c, for the first time. We present these 
interferometric observations along with additional spectroscopic radial velocity 
measurements of the components. We performed a simultaneous orbit fit to the 
interferometric visibilities and the radial velocities of Haro l-14c. Based on a 
statistical analysis of the possible orbital solutions that fit the data, we deter- 
mined component masses of M x = 0.96± f 8 M Q and M 2 = 0.33^02 M & for the 
primary and secondary, respectively, and a distance to the system of 111 pc. 
The distance measurement is consistent with the close distance estimates of the 
Ophiuchus molecular cloud. Comparing our results with evolutionary tracks 
suggests an age of 3-4 Myr for Haro l-14c. With additional interferometric mea- 
surements to improve the uncertainties in the masses and distance, we expect 
the low-mass secondary to provide important empirical data for calibrating the 
theoretical evolutionary tracks for pre— main-sequence stars. 

Subject headings: binaries: spectroscopic — binaries: visual — stars: fundamen- 
tal parameters — stars: individual (Haro l-14c) — stars: pre— main-sequence - 
techniques: interferometric 



1. Introduction 

For masses smaller than ~ 1M Q , uncertainties in the pre— main-sequence (PMS) evolu- 
tionary tracks are sufficiently great that the usual method of estimating the mass and age of 
a young star, by its location in the H-R diagram, can produce a scatter as large as a factor of 
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2 in mass and 3 in age depending on the tracks used (e.g.. lSimonll200ll ; lMathieu et al.ll2007l ). 
As a result, the mass spectrum of stars produced in a star- forming region, the distribution of 
masses in binaries, and the region's star-forming history are imprecisely known. Measuring 
reliable masses, luminosities, and effective temperatures of young stars contributes empirical 
data for testing calculations of young star evolution. Currently, there are few PMS stars 
with reliably mea sured masses; in particular, there are none below 0.4 M & with a precision 
better than 10% Jffillenbrand k Whitelliooi iMathieu et alil2007h . 



Mapping the orbital motion of a binary system provides a way to measure the dynamical 
masses of PMS stars. Observations of a system as both a visual and spectroscopic binary 
yield the component stellar masses as well as the distance to the system. Spectroscopic 
radial velocity measurements of PMS binaries are typically limited to systems which have 
periods smaller than a few years. At the distance of nearby star forming regions like Taurus 
and Ophiuchus, this range of orbital periods corresponds to angular separations smaller than 
a few tens of milli-arcseconds (mas). As a result, there are few PMS spectroscopic binaries 
that can be resolved spatially using the current generation of high resolution instrumentation. 
Steffen et all fhoOll ) applied this approach to the PMS binary NTTS 045251+3016 by com- 
bining single- and double-lined radial velocity measurements with spatially reso lved observa- 



tions obtained with the Fine Guidance Sensors on the Hubble Space Telescope. iBoden et al. 



( 120051 ) used the Keck Interferometer and published spectroscopic observations to determine 
the masses of the PMS binary HD 98800B. Precise mass measurements of even a few systems 
provide critical data for calibrating the evolutionary tracks. Additionally, the distances de- 
termined from the joint orbital solutions provide an accurate location of the system within 
the relatively spread-out star forming regions. 



Haro l-14c (HBC 644) is a PMS binary in the Ophiuchus star forming region. iReipurth et al 



( 120021 ) first identified the system as a single-lined spectroscopic binary through observations 
in the visible part of the spectrum. By observing in the infrared, where the light ratio be- 
tween the primary and the cooler secondary is more favorable, Simon & Prato (2004, SP04) 
measured the radial velocity amplitude of the secondary and determined a secondary to pri- 
mary mass ratio of 0.310 ±0.014. In this paper we present spatially resolved observations of 
Haro l-14c obtained with the Keck Interferometer. We also present additional spectroscopic 
measurements to improve the precision of the velocity amplitude of the secondary. These 
observations are described in § |2j In § [3j we present results from a simultaneous orbit fit 
to the interferometric visibilities and the spectroscopic radial velocities. In § HI we discuss 
the derived distance to Haro l-14c and compare the dynamical masses with predictions from 
evolutionary tracks. A summary of the results appears in § [51 
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2. Observations of Haro l-14c 
2.1. Interferometric Observations 

The Keck Interferometer combines the light of t he two 10-meter Keck telescopes and has 



a bas eline of 85-meters, oriented 38° east of north (IColavita et al.l 120041 ; IWizinowich et al. 



20041 ) . The observable for the interferometer is the visibility or fringe contrast of the source. 
For a binary star where the angular diameter of each component is unresolved, the normalized 
visibility V 2 at a wavelength A is 

V 2 _ 1 + P + V cos[27r/A(nAa + vA5)} 

where (u,v) are the baseline components projected on the sky, (Aa,A<5) are the binary 
separations in RA and DEC, and / is the component flux ratio (e.g. Boden 1999). At the 
distance of the Ophiuchus star forming region, we expect the angular diameters of the stellar 
components of Haro l-14c to be unresolved by the Keck Interferometer. 

We used the Keck Interferometer in its V 2 science mode, operating in the K-band. We 
observed Haro l-14c on five nights between June 2004 and May 2007. To calibrate the 
fringe visibilities, we interspersed observations of unresolved single stars between those of 
Haro l-14c. An unresolved point source will have a visibility amplitude of 1.0. We selected 
the calibrators using the getCal planning tool distributed by the Michelson Science Center 
(MSC). The properties of Haro l-14c and the selected calibrators are given in Table [T] which 
lists the RA and DEC coordinates, spectral type, V and K magnitudes, parallax, separation 
on the sky from Haro l-14c, estimated angular diameter, and epochs of observation. The 
angular diameters were calculated using the fbol routine in the getCal package which fits 
a blackbody model to spectrophotometric data. Input to the routine included B and V 
magnitudes from the All-sky Compiled Catalog of 2.5 million stars (Kharchenko 2001), J, 
H, K band magnitudes from the Two Micron All-Sky Survey (2MASS) Catalog (Cutri et al. 
2003), and parallaxes from the Hipparcos Catalog (Perryman et al. 1997). The last three 
calibrators listed in Table [T] were observed as part of other programs but were used in this 
paper to improve the time sampling of calibrators surrounding the Haro l-14c observations. 

We calibrated the raw visibilities using the wbcalib package developed by the MSC. In 
doing so, we averaged the raw V 2 data over scan lengths of 300-500 sec. We used the standard 
deviation of the raw measurements as an estimate of the uncertainties. The wbcalib package 
computes the system visibility as a function of time based on the unresolved calibrator 
observations and applies this correction to the science data. We also applie d a correction 



to ac count for a flux-dependent bias between the target and the calibrators (lAkeson et al. 



20071 ). Table[2]lists the Modified Julian Date (MJD), UT date and time, observed wavelength, 
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calibrated V 2 and error, and the u and w projected baselines during the times of observations. 
The calibrated visibilities of Haro l-14c measured by the Keck Interferometer in 2004-2007 
are plotted in Figure [TJ Figure [2] shows an example of the wf-plane coverage obtained on 
2004 June 2. 



2.2. Velocity Measurements 



Near-infr ared spectra were ob t ained at the Keck II telescope using the NIRSPEC facility 
spectrograph JMcLean et al.lllQQsl . I2OO0I ) . NIRSPEC employs a 1024x1024 ALADDIN InSb 
array detector with a plate scale of 0.14" in the dispersion direction and 0.19" in the cross- 
dispersion direction. The observations were centered at 1.555 fim in the H-band. A 0.288" 
(2 pixel) by 24" slit yielded high resolution (R=30,000) spectra in several cross-dispersed 
orders. For the new observations, 2—4 minute exposure times in nodded ABBA sets were 
typical; the seeing was 0.6—0.8". Further information on the previously published NIRSPEC 
data can be found in SP04. 

We used only the central order, 49, in our analysis because it covers a region rich in 
both atomic and molecular lines suitable for the characterization of a range of spectral type 
components. Furthermore, while this spectral region lacks contamination from strong ter- 
restrial absorption lines, more than 10 relatively st rong OH night sky em ission lines span 
the order, providing suitable dispersion calibration (IRousselot et al.ll2000l ). All data reduc- 
tion was performed using the REDSPEC packag^j], designed specifically for the reduction of 
NIRSPEC data. 

Table 3 lists the UT date, Modified Julian Dates (MJD), velocities and uncertainties 
for the primary and secondary, respectively, and the orbital phase determined from the sub- 
sequent orbital fit (§3.) The first 6 lines in Table 3 represent a reanalysis of the observations 
reported in SP04. The last 2 lines refer to the new observations made after the publication 
of SP04. 

To measure the spectral types and velocities of the components we analyzed all the 
spectra, including the spectra used in SP04, as follows. We first used TODCOR, the two- 
dimensional r _cross I con^ljtion algori thm written for analysis of double-lined spectroscopic 
binaries by 



( jPrato et al. 



Zucker fc Mazehl (119941 ) and our suite of observed stellar spectral templates 
20021 ) to identify the spectral templates providing the best match in spectral 



type and rotational broadening. HR 8085, a K5 V, and GL 15A, a Ml. 5 V, both rotationally 



http: / /www2. keck. hawaii.edu/inst/nirspec/redspec/index. html 
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broadened to 12 km s l , provided the best matches for the primary and secondary of Haro 
l-14c, consistent with the result reported in SP04. 

We chose not to use the observed stellar templates to measure the component velocities, 
as we had in SP04, because to do so would have propagated velocity uncertainties associated 
with the templates to Haro l-14c. Ins tead, we used a suite of theoretical spectra calculated 



from updates of the NextGen models (lHauschildt et al.lll999l ). These calculations represent 



main sequence stars with solar abundance in the mass range 0.1 to 1.6 M & . We selected 
theoretical spectra that provided excellent matches to the observed spectra of HR 8085 and 
GL 15A. We then used these two corresponding spectra, with a velocity correction to the 
heliocentric frame appropriate for each target observation, to measure the radial velocities of 
the primary and secondary components of Haro l-14c. We applied Wilson's (1941) approach 
to our 8 pairs of primary and secondary velocities to determine the center-of-mass velocity 
and its uncertainty, 7 = —9.38 ± 0.39 km s _1 . This value of 7 and the component velocities 
we measured from the infrared spectra are i n the laboratory frame of the theoretical spectra, 



which is not necessarily the frame in which iReipurth et al.l (120021 ) measured the velocity of 



the Haro l-14c primary. Their 56 measurements yield a very accurate value of the center- 
of-mass velocity of 7 = —8.71 ± 0.07 km s _1 . We chose therefore to shift our velocity 
measurements by the difference between our and Reipurth et al.'s value of 7, —0.67 ± 0.39 
km s -1 , effectively transforming them into the frame of Reipurth et al. These values, and 
their propagated uncertainties, are reported in Table 3. The velocity uncertainties of the 
secondary, which are particularly important in the determination of the secondary/primary 
mass ratio, are about half the values reported in SP04. 



3. Orbital Fit to Haro l-14c 

Velocity measurements of a double-lined spectroscopic binary yield the orbital period, 
eccentricity, time of periastron passage, longitude of periastron, and the velocity semi- 
amplitudes of the primary and secondary (P, T, e, lu, Ki, and K 2 ). A visual binary yields 
P, T, e, u, a("), i, and Q, where the last three parameters are the semi-major axis in arc- 
seconds, the inclination, and the position angle of the line of nodes. Combining the visual 
and spectroscopic orbital parameters gives the masses of the primary and secondary, 

. . _ 1.036 x lO- 7 ^ + K 2 fK 2 P(\ - e 2 f' 2 

M x — — [Z 

strict 

A/f _ 1-036 x lO- 7 ^ + K 2 fK x P{\ - e 2 ) 3 / 2 
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where K\ and K% are in km s -1 , P in days, and Mi and M 2 in M & (jFernid |2000| ). The 
distance to the binary is determined by comparing the angular size of the visual orbit with 
the physical scale of the spectroscopic orbit, 



d 



9.198 x KT 5 ^! + K 2 )P(l - e 
a(") sini 



au/2 



(4) 



where the distance d is in pc. 



3.1. Formal Orbital Analysis 

We performed a simultaneous orbit fit to the i nterferometer v is ibiliti es and spectroscopic 



radial velocities presented in this paper and in iReipurth et al.l (120021 ). The orbit fitting 
procedure starts with initial estimates of the orbital parameters (P, T, e, u, a, i, Q, K\, 
K 2 , 7) and the K-band flux ratio and follows a Newton-Raphson approach to adjust these 
parameters to their best-fit values. During each iteration of the routine, the values of (Aa, 
AS) are computed at the times of the interferometer observations from the given set of 
orbital parameters. The program then calculates model visibilities according to Equation [TJ 
assuming that the angular diameters of the components are unresolved. The model radial 
velocities are determined from 

V\ = +K 1 [e cos oj + cos (u + u)\ + 7 (5) 

v 2 = — K 2 [e cos oj + cos (v + u)] + 7 (6) 

where the true anomaly u is computed at the times of the spectroscopic observations from 
the given set of orbital parameters. The program adjusts the initial parameters to minimize 
the x 2 between the measurements and model values, which is given by 

2 _ (Kncas ~ ^fit) 2 + (^lmeas ~ ^lfit) 2 + (^2mcas ~ ~^2fit) 2 

where the 'meas' and 'fit' subscripts refer to the measured values and the model fit, respec- 
tively. 

Table H] lists the orbital parameters, component masses, and distance determined from 
the simultaneous orbit fit. The last two rows show the \ 2 °f the fit and the reduced xl f° r 
83 degrees of freedom. The best-fit orbit is overplotted in the V 2 curves shown in Figure [TJ 
Figure [3] shows the orientation of the binary on the plane of the sky. The epochs of the Keck 
Interferometer observations are marked along the orbit. The radial velocity measurements 
and best-fit velocity curves are shown in Figure HI 
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The spectroscopic orbital parameters (P, T, e, lu, Ki, K2, 7) of Haro l-14c are well- 
determined by the single- and double-lined radial velocities. However, the limited sampling 
of the visual orbit provided by the Keck Interferometer observations produces a large uncer- 
tainty in the inclination. We examined the distributions of a and i allowed by the current 
set of V 2 measurements by performing a two-dimensional Monte Carlo search. During each 
iteration of the search, we selected values of a and i at random and minimized the remain- 
ing orbital parameters through the standard Newton-Raphson approach. We accumulated 
a total of 10,000 possible solutions within the 3-<r confidence interval (A% 2 = 9 from the 
minimum \ 2 value). The resulting x 2 distribution is shown in Figure [51 The long tail of so- 
lutions that extends to small inclinations produces extremely large component masses when 
combined with the spectroscopic parameters. The uncertainties for a, i, Q, Mi, M2, and d 
quoted in Table H] represent the bounds of the 1-a (Ax 2 = 1) confidence regions. 

The results from the Monte Carlo search also allow us to examine the uncertainties in 
the binary separation derived from the epochs of the Keck Interferometer observations. For 
each of the 10,000 orbital solutions found, we computed the binary separations at the times of 
the interferometer observations and overplotted the 1-a confidence intervals in the orbit plot 
in Figure [3J The elongated error ellipses arise from the limited wf-coverage sampled during 
the interferometer observations (see Figure [2]). Given the low altitude of Haro l-14c at the 
latitude of the Keck Observatory, the short 1-3 hour observing sessions do not provide much 
rotation of the projected interferometer baseline on the sky. Essentially, the narrowness of 
the error ellipses shows the exceptional resolution along the direction of the interferometer 
baseline. Because we fit all of the visibilities from each epoch simultaneously with the 
spectroscopic radial velocities, additional measurements with the Keck Interferometer at 
different locations along the orbit will help reduce the range of allowed orbital solutions and 
hence reduce the size of the error ellipses in the binary separation for every epoch (see § 13.31 
for more details). 

3.2. Masses and Distance from a Statistical Analysis of Orbital Solutions 

To investigate the distribution of masses and distance allowed by the interferometric and 
spectroscopic data, we performed a more exhaustive 6-dimensional Monte Carlo search. We 
randomly selected values for all of the orbital parameters on which the masses and distance 
depend (P, e, a, i, Ki, and K 2 ) while minimizing the fit to the remaining parameters through 
the standard Newton-Raphson approach. We accumulated a total of 10,000 possible solutions 
within the A% 2 = 9 confidence interval. 

The top row of Figure M shows cross-cuts through the x 2 surface for the component 
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masse s and distance versus the inclination. Following the approach outlined in lSchaefer et al. 



( 120061 ). we constructed histograms of the allowed masses and distance where each solution 
is weighted by its \ 2 probability. As seen in these histograms, shown in the bottom row of 
Figure El the mass and distance distributions are strongly peaked. We present the median 
values of these distributions in Table [5j The range of the quoted uncertainties include 34% 
of the values on each side of the median. In computing these values, we applied an arbitrary 
upper mass cut-off of 2.0 M Q for the primary, since masses larger than this value are clearly 
not representative of the spectral type and stellar luminosity of Haro l-14c. Compared to 
the formal x 2 intervals in Table HI the statistical analysis substantially reduces the size of 
the error bars. It is worth noting that because of the sharp drop in the mass distributions 
on the low-mass sides, improved mass values are not likely to be much smaller than the 
present values, but may be larger. In the following sections we refer to the median masses 
and distance estimates given in Table as we believe these values more accurately represent 
the asymmetric parameter distributions. 



3.3. Comments on the Reliability of the Orbital Fit and Expectations for the 

Future 

Obtaining interferometric observations that are well distributed across the orbit is crit- 
ical for determining a reliable measurement of the inclination. Since the inclination is de- 
termined from the offset of the primary star from the focal point of the apparent orbit, 
measuring an accurate value requires good orbital sampling along both the major and minor 
axes of the apparent ellipse. For the orbit of Haro l-14c, we have fairly good coverage along 
the major axis of the ellipse, but we have only one marginal measurement along the direction 
of the minor axis. During the observations on 2005 April 20, we experienced unusually low 
instrument throughput on Haro l-14c which made it difficult to lock onto the fringes and 
consequently degraded the quality of the V 2 measurements. Because these marginal mea- 
surements sample a critical location along the orbit, a precise value of the inclination has 
not yet been determined. 

To investigate what we can expect from future observations with the Keck Interfer- 
ometer, we simulated V 2 measurements for 2008 June 1 and 2012 May 1. These dates 
were selected to provide good sampling of the Haro l-14c orbit at times when the binary is 
observable with the Keck Interferometer. To simulate the interferometer measurements, we 
computed the expected V 2 based on the predicted binary position at 20 minute intervals over 
2.3 hours centered on transit. We added Gaussian noise to the simulated V 2 and recomputed 
the simultaneous orbit fit based on the measured and simulated V 2 and spectroscopic radial 
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velocities. As seen in Figure 0, the additional measurements would significantly reduce the 
uncertainties in the orbit fit. Moreover, the simulated data yield formal \ 2 uncertainties 
in the masses and distance that are smaller by more than a factor of 2 compared with the 
current values in Table HI 



Discussion 



4.1. Distance to Haro l-14c 



The orbital measurements yield a dynamical estimate of the distance to Haro l-14c of 
d = llliig P c - This result is consiste nt with the close dista nce estimates of th e Ophi uchus 



dark clou ds (~ 120 pc) determ ined by iKnude fc Hpgl (119981 ) and Ide Geus et all (119891 ) . Ad- 
ditionally, lLoinard et al.l (120081 ) find a similar mean distance of 120 ± 4 pc to the Ophiuchus 
core using the Very Long Baseline Array (VLBA) to measure the parallaxes of radio-emitting 
PMS stars. The distance of ~ 116 pc determined from the formal orbital parameters (Ta- 
ble HJ) is slightly closer to the VLBA value than our statistical estimate, although both are in 
agreement within the l-er uncertainties. Ultimately, accurate distances that map the spatial 
extent of the Ophiuchus clouds will be necessary for determining the luminosities, and hence 
ages, of young stars in the region. 



4.2. Absolute Magnitudes and Effective Temperatures 

In addition to the distance estimate, the interferometric measurements yield the flux 
ratio of the components in the K-band. Because we have insufficient information to determine 
the luminosities, in § 14.31 we compare our results with the theoretical PMS evolutionary 
tracks using H-R diagrams plotting the absolute magnitude at K, Mk, versus the effective 
temperature, T e g. The primary and secondary absolute magnitudes are given by 

M K1 or 2 = m K1 or 2 - 5 log d + 5 - A K (8) 

where d is the distance in pc, lll^U pc, and Ax the extinction at K. To calculate the 
extinction at V we use 



A v = 13.83(J-#) 



observed 



8.29(H — K) 



observed 



- 7.43 



(9) 



derived by iPrato et al.l (120031 ) from the reddening relations used b y iMeyer et al.l (Il997l ). 
Then, with A x = (0.60/A) L75 A y for 0.9 < A < 6/im jTokunagalboOoh . A K = 0.10A V . 
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The apparent magnitudes of the Haro l-14c binary available in the 2MASS Catalog at 
J, H, and K are 8.86±0.03, 8.00±0.05, and 7.78±0.03 mag, respectively. With the measured 
binary flux ratio at K (Table 4), the apparent magnitudes of the primary and secondary are 
m Ki = 8.04±0.04 and rriKi = 9.47±0.14 mag. Using the 2MASS values (which apply to the 
light of both components) in equation^ we obtain a preliminary estimate of Ay ~ 2.64±0.96 
mag. We used the PMS models calculated by Siess et al. (2000) to estimate the effect of 
the two components and derive a corrected visual extinction of Ay ~ 2.9 ± 1.0 mag. Using 
the relation Ak = O.lOAy, we find A^ = 0.29 ± 0.10 mag. The absolute magnitudes of 
the primary and secondary are therefore Mki = 2.52 ± 0.04 and Mk2 = 3.95 ± 0.14 mag. 
The uncertainties include only the errors in the measured K-band flux ratio and the 2MASS 
photometry but not the uncertainties in either the distance modulus or extinction because 
their values apply equally to both components. We will display the effects of the uncertainties 
in the distance modulus and extinction on the H-R diagrams presented in §4.31 

Based on the observed spectral templates, the best fitting spectral types for the pri- 
mary and secondary of Haro l-14c are K5 and Ml. 5, respectively. To estimate the effec- 
tive temperatures of the pri mary and secondary, we used the temperature scale adopted by 



Hillenbrand fc White! (120041 ) and derive values of 4395 ±20 K for the primary and 3548 ±245 



K for the secondary. We assume uncertainties of ± 1 spectral subclass. 



4.3. Comparisons with the Evolutionary Tracks 

The dynamical masses of M 1 = 0.96^qo§ M and M 2 = 0.33^oo2 M Q are in the range 
where precise values will make a significant contribution to testing the theoretical tracks of 
pre-main sequence stellar evolution. This is particularly the case for the low mass secondary. 
Moreover, the large difference in the mass of the primary and secondary will facilitate ap- 
plying the coevality test to the tracks and estimating the age of the system. Based on our 
statistical analysis of possible orbits, the precision of the present values of the primary and 
secondary masses is about 28% and the precision in the distance is 17%. In the following 
sections we compare our mass estimates with the evolutionary tracks of Baraffe et al. (1998, 
2002; BCAH) and Siess et al. (2000; SDF). The evolutionary tracks are plotted in Figure [HJ 
where we highlight the tracks interpolated to our dynamical mass estimates of the primary 
and secondary of Haro l-14c. We overplotted the location of the effective temperatures and 
absolute K-band magnitudes determined for the primary and secondary ( §4.21) . We discuss 
the comparisons between each set of models in more detail below. 
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4.3.1. SDF Models 



The SDF tracks are available at subsolar, solar, and supersolar metallicities (Z=0.1, 0.2, 
0.3, respectively). We retrieved the absolute K-b and magnitudes from the S DF database 
that were computed using the temperature scale of iKenyon fc Hartmannl (119951 ). The tracks 
suggest that the Haro l-14c binary is 3-4 Myr old. Arguably, our measurements fit the tracks 
for each of the three abundances within the uncertainties, but with a slight preference for 
the solar or supersolar tracks. 



4.3.2. BC AH Models 

The BCAH tracks are available for solar and subsolar ([m/H]=-0.5) metallicities. For 
masses below 0.7 M , we used the tracks calculated with a mixing length parameter of 1.0; 
for greater masses we used the tracks with a mixing length of 1.9. The solar and subsolar 
tracks indicate an age for the binary of ~ 3 Myr, or slightly older. Agreement with the mass 
and coevality appears to be slightly better on the subsolar tracks. 



4-3.3. Summary of Comparisons 



Across the range of metallicities considered, both the tracks of BCAH and SDF indicate 
an age of 3-4 Myr for the Haro l-14c binary. Sliding the components up and down according 
to the distance and extinction uncertainties in the H-R diagrams would change the binary 
age but not the qualitative agreement with the tracks. Overall, the agreement between the 
evolutionary tracks and our measurements seems to be better for the models of SDF than 
for BCAH, especially when considering the coevality criterion. 

The evolutionary tracks of SDF and BCAH allow for the possibility of solar, supersolar, 
and subsolar metallicities for Haro l-14c. However, based on a comparison with the solar 
metallicity spectral template, HR 8085, we argue that a subsolar or supersolar metallicity 
for Haro l-14c is unlikely given the strength of the Fe lines acr oss order 49 of the NIR- 
SPEC spectrum (see Figure [9]). Additionally, ISantos et al.l (120081 ) find a metallicity that is 
indistinguishable from solar for the Ophiuchus star-forming region. 

Our age estima te of 3-4 Myr for Har o l-14c is slightly older than the median age of 2.1 
Myr determin e d by IWilking et al.l (120051 ) for the p Ophiuchus molecular cloud. However, 
Wilking et al.l (120051 ) use a distance of 150 pc in determining the stellar luminosities in their 
sample. If a large portion of the stars in Ophiuchus are located at the closer distance estimate 
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of 120 pc, the net effect would be to decrease the intrinsic luminosities and increase the age 
estimate of stars in the region. 

From our preliminary analysis, it appears that a major limiting factor in making com- 
parisons with the PMS evolutionary tracks lies in our ability to estimate the effective tem- 
peratures of the component stars. Based on the spectral type s of the best-fitting observed 
spectral templates, we used the temperature scale adopted by Hillenbrand fc White! ( 120041 ) 
to estimate T e g for the components of Haro l-14c. Temperatures of PMS stars determined in 
this way are not only affected by the ± 1 spectral subclass uncertainties from the fit, but are 
also subject t o differences in the adopted t emperature sca l e. Fo r inst ance, the temperature 
scales used by iKenyon &: Hartmannl (119951 ) , iLuhman et al.l (120031 ) , and iHillenbrand &: White 
(120041 ) vary by more than 200 K for stars of spectral type M2 and later. Possible improve- 
ments in the T e g measurements of PMS stars could b e obtained from fitting sy nthetic spectral 
templates to high resolution near-infrared spectra (IDoppmann fc Jaffd 120031 ) . however, this 
may prove difficult for binaries where the spectra of the rotationally broadened components 
are blended. 



5. Conclusions 

1. Using the Keck Interferometer, we spatially resolved the orbit of the pre-main se- 
quence binary, Haro l-14c. We also added two new spectroscopic measurements of the radial 
velocities of the primary and secondary. 

2. We computed a simultaneous orbit fit to the interferometer visibilities and the spec- 
troscopic radial velocities. Based on a statistical analysis of the orbital parameters that fit the 
data, we determined component masses of Mi = 0.96^0^ M Q and M 2 = 0.331q;q2 M q . The 
secondary in Haro l-14c is one of the lowest mass PMS stars with a dynamical measurement 
of its mass. 

3. Our orbital measurements also yield a dynamical distance to Haro l-14c of lllljj pc. 
This result is consistent with the close distance estimates to the Ophiuchus molecular cloud. 

4. Comparing our results with the evolutionary tracks of BCAH and SDF suggests an 
age of 3-4 Myr for Haro l-14c. 

5. The tracks of SDF appear to be more consistent with our measurements of Haro 
l-14c than those of BCAH, especially when considering the coevality of the components. 

6. Additional interferometric observations of Haro l-14c at well-chosen times during the 
orbit will substantially improve the estimate of the orbital inclination. Not only will this 
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provide precise values for the dynamical masses of the components, but it will also improve 
the distance measurement for the system. An accurate distance is needed to derive the 
luminosities and place the stars on the H-R diagram. The remaining issue that needs to be 
addressed in comparing dynamical mass measurements with theoretical evolutionary tracks 
lies in determining accurate effective temperatures for PMS binaries. 
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Table 1. Calibrators for Keck Interferometer Observation of Haro l-14c 



Object 


Coordinates 


SpT 


ya 


K b 




Sep 


Ang Diam 


Epoch d 














(mas) 


(°) 


(mas) 




Haro l-14c p 


16:31:04.36 


-24:04:32.7 


K5V 


12.3 


8.0 








1,2,4,5 


Haro l-14c s 


16:31:04.36 


-24:04:32.7 


M1.5V 




9.5 








1,2,4,5 


HD 143096 


15 59 25.863 


-30 19 57.046 


K0V 


9.6 


7.4 


19.2 


9.1 


0.17±0.07 


1 


HD 143955 


16 04 00.746 


-19 06 15.841 


G5V 


9.2 


7.5 


15.3 


8.1 


0.14±0.04 


3,4 


HD 147935 


16 25 49.189 


-27 49 09.193 


G5V 


9.2 


7.6 


17.5 


2.6 


0.14±0.03 


1,2,4,5 


HD 151692 


16 49 53.157 


-24 26 48.976 


K3V 


9.6 


6.9 


29.1 


4.3 


0.25±0.12 


1,4,5 


HD 152400 


16 54 14.962 


-28 06 53.799 


G6V 


9.2 


7.7 


12.9 


6.5 


0.13±0.02 


1 


HD 155006 


17 09 50.413 


-21 47 31.502 


G6V 


9.2 


6.4 


10.1 


9.7 


0.31±0.11 


5 


HD 147284 


16 21 55.439 


-24 59 28.690 


G3V 


8.8 


7.3 


13.5 


2.3 


0.16±0.05 


3 


HD 148376 


16 28 02.566 


-13 06 10.488 


F7V 


8.7 


7.2 


9.8 


11.1 


0.15±0.02 


2 


HD 153022 


16 57 27.164 


-12 14 04.254 


G8V 


8.8 


5.8 


4.6 


13.7 


0.43±0.20 


2 



a V-band magnitudes from the All-sky Compiled Catalog of 2.5 million stars (Kharchenko 2001) 
b K-band magnitudes from the 2MASS Catalog (Cutri et al. 2003) 
c Parallaxes from the Hipparcos and Tycho Catalog (Perryman et al. 1997) 

d Epoch: 1 = 2004 June 2; 2 = 2005 April 20; 3 = 2006 April 4; 4 = 2006 May 17; 5 = 2007 May 1 
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Table 2. Haro l-14c V 2 Measurements 



MJD 


Date 


U 1 


A 


V 2 


ay 2 


u 


V 










(/mi) 










(m) 


(mj 


53158 


,37862 


2004 Jun 2 


09:05:12 


2. 


18 


0. 


.795 


0. 


.152 


55. 


.98 


54.70 


53158. 


,40510 


2004 Jun 2 


09:43:20 


2. 


18 


0. 


.873 


0. 


.123 


53 


.89 


50.95 


53158. 


.41680 


2004 Jun 2 


10:00:11 


2. 


18 


0. 


.881 


0. 


090 


52. 


.49 


49.35 


53158. 


.42913 


2004 Jun 2 


10:17:56 


2. 


18 


0. 


.947 


0. 


.083 


50. 


.70 


47.72 


53158. 


.44104 


2004 Jun 2 


10:35:06 


2. 


18 


0. 


.776 


0. 


.119 


48. 


.68 


46.19 


53158. 


.45318 


2004 Jun 2 


10:52:34 


2. 


18 


0. 


.686 


0. 


.089 


46. 


.34 


44.71 


53158. 


,46304 


2004 Jun 2 


11:06:46 


2. 


18 


0. 


.645 


0. 


.070 


44. 


.24 


43.56 


DOlDO 


,4 ( o4U 


2004 Jun 2 


11:21:41 


2. 


18 


U. 


4DZ 


u 


U (v 


41. 


.85 


42.41 


53480. 


,57093 


2005 Apr 20 


13:42:08 


2. 


.18 


0. 


.295 


0. 


.162 


46. 


.41 


44.75 


53480. 


,61274 


2005 Apr 20 


14:42:20 


2. 


.18 


0. 


.252 


0. 


223 


36. 


.41 


40.27 


53839 


,58124 


2006 Apr 14 


13:56:58 


2. 


.15 


0. 


.824 


0. 


.055 


47. 


.75 


45.57 


53839 


.60523 


2006 Apr 14 


14:31:32 


2. 


.15 


0. 


.946 


0. 


.067 


42. 


.64 


42.77 


53839 


.62265 


2006 Apr 14 


14:56:37 


2. 


.15 


0. 


.719 


0. 


.028 


38. 


32 


40.96 


53839 


.63442 


2006 Apr 14 


15:13:33 


2. 


.15 


0. 


.495 


0. 


.038 


35. 


.14 


39.85 


53872. 


.45111 


2006 May 17 


10:49:35 


2. 


.15 


0. 


.498 


0. 


.037 


53. 


.79 


50.82 


53872. 


,46430 


2006 May 17 


11:08:35 


2. 


.15 


0. 


.646 


0. 


.042 


52. 


16 


49.03 


54221. 


,47782 


2007 May 1 


11:28:03 


2. 


18 


0. 


.890 


0. 


.098 


55. 


.40 


53.31 


54221. 


.48940 


2007 May 1 


11:44:43 


2. 


18 


0. 


.983 


0. 


.035 


54. 


.42 


51.67 


54221. 


,51506 


2007 May 1 


12:21:40 


2. 


18 


0. 


991 


0. 


.045 


51. 


.25 


48.18 


54221. 


,54398 


2007 May 1 


13:03:19 


2. 


18 


0. 


.898 


0. 


.019 


46. 


.07 


44.55 


54221. 


.55555 


2007 May 1 


13:19:59 


2. 


18 


0. 


.818 


0. 


019 


43. 


.56 


43.21 


54221. 


.56694 


2007 May 1 


13:36:23 


2. 


18 


0. 


.769 


0. 


.035 


40. 


.86 


41.97 
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Table 3. Haro l-14c Radial Velocity Measurements 



UT Date 


MJD 


Vl 




v 2 


^2 


Phase 






(km s _1 ) 


(km s -1 ) 


(km s _1 ) 


(km s _1 ) 




2001 Jun 1 


52061.51 


-4.34 


0.59 


-17.51 


1.16 


0.308 


2002 Jul 17 


52472.28 


-18.13 


0.60 


15.00 


1.16 


0.002 


2003 Feb 8 


52678.67 


-6.14 


0.62 


-17.49 


1.07 


0.350 


2003 Aug 10 


52861.32 


-9.68 


0.54 


-6.00 


1.14 


0.659 


2003 Sep 8 


52890.20 


-10.59 


0.57 


-4.27 


0.94 


0.707 


2004 Jan 26 


53030.69 


-21.67 


0.61 


23.74 


1.08 


0.945 


2007 Apr 30 


54220.48 


-21.17 


0.58 


26.86 


1.07 


0.954 


2007 Aug 9 


54321.26 


-3.55 


0.56 


-22.85 


1.36 


0.124 



Table 4. Haro l-14c Orbital Parameters from a Formal Analysis 



Parameter Value 

P (days) 592.11 ± 0.12 

T (MJD) 53655.51 ± 0.86 

e 0.6211 ± 0.0070 

u(°) 232.4 ± 1.2 

a (mas) 13.0j£? 

75.9^3 

m 66+H 

K x (km/s) 8.73 ± 0.13 

K 2 (km/s) 25.53 ± 0.56 

7 (km/s) -8.788 ± 0.076 

/(k) o.268ig;82g 

M 1 (M Q ) 0.975±8:8K 

M 2 (M ) 0.333^| 



-38 
-21 



d (pc) 116. 
X 2 101.86 
Xl 1-227 
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Table 5. Best-fit Masses and Distance from a Statistical Analysis of Orbital Solutions 



Parameter 


Value 


M X {M & ) 
M 2 (M ) 
d (pc) 


Qfil +0-265 

111 + 19 
lll._ 18 
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Fig. 1. — Calibrated visibilities of Haro l-14c measured with the Keck Interferometer on UT 
2004 June 2, 2005 April 20, 2006 April 14, 2006 May 17, and 2007 May I. Overplotted in each 
panel is the best-fit binary orbit computed from a simultaneous fit to the V 2 measurements 
and the spectroscopic radial velocities. 
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Fig. 2. — uv -coverage sampled by the Keck Interferometer observations of Haro l-14c on UT 
2004 June 2. 
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Fig. 3. — Best-fit orbit for Haro l-14c. The large squares indicate the orbital position at the 
epochs of the Keck Interferometer observations. The elongated error ellipses represent the 
1-er confidence intervals (A% 2 = 1) determined from the Monte Carlo search of the orbital 
parameter space. The narrowness of the ellipses (0.1-0.3 mas) reflects the precision along 
the direction of the interferometer baseline. 
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Phase 

Fig. 4. — Radial velocity as a function of phase for Haro l-14c. The crosses represent the 
Reipurth et al. (2002) measurements of the primary velocity. The open symbols and la 
uncertainties indicate the radial velocities for the primary (circles) and secondary (squares) 
measured using NIRSPEC. The solid line shows the best fit to the primary star data, and 
the dotted line for the secondary. The horizontal dashed line indicates the center-of-mass 
velocity of the system. The abscissa spans a 20% redundancy in phase for clarity. 
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100 




a (mas) 



Fig. 5. — x 2 surfaces determined through the two-dimensional Monte Carlo search for orbital 
solutions that fit the current data set. The color codes correspond to the 1 o (red), 2 o 
(green), and 3 a (blue) confidence intervals (A% 2 = 1, 4, and 9, respectively). 
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Fig. 6. — Top row: Crosscuts through the \ 2 surface for M\ (left), M 2 (middle), and distance 
(right) versus the inclination based on the results from the 6-dimensional Monte Carlo search 
described in § 13.21 The color codes correspond to the 1 a (red), 2 a (green), and 3 er (blue) 
uncertainty intervals (A% 2 = 1, 4, and 9, respectively). We added an additional 1,000 
solutions in each of the 1 a and 2 o intervals to better define these regions visually in the 
plots; these added solutions were not used in the statistical analysis of the distributions. 
Bottom row: Weighted distributions for Mi (left), M 2 (middle), and distance (right). In 
constructing the histograms, each orbital solution was weighted by its x 2 probability. The 
dashed line indicates the solutions removed by the 2 M & upper mass cut-off for the primary. 
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Fig. 7. — Modeled orbit fit based on simulated Keck Interferometer data for 2008 June 1 and 
2012 May 1. These dates were selected to provide good sampling of the Haro l-14c orbit at 
times when the binary is observable with the Keck Interferometer. The expected improve- 
ment in the orbital fit with additional measurements is clearly apparent when comparing to 
the current fit in Figure [3j 




Fig. 8. — Top row: Evolutionary tracks computed by SDF. The tracks (dashed lines) cover 
the mass range from 0.1 to 1.6 M , plotted at 0.1 M intervals. The isochrones (solid lines) 
are plotted at ages of 1, 3, and 10 Myr. The three panels show the SDF tracks computed at 
subsolar (Z=0.01; left), solar (Z=0.02; middle), and supersolar (Z=0.03; right) metallicities. 
Bottom row: Evolutionary tracks computed by BCAH. The tracks (dashed lines) cover the 
mass range from 0.1 to 1.4 M , plotted at steps of 0.1 M out to 1.0 M and steps of 0.2 M 
out to 1.4 M & . The isochrones (solid lines) are plotted at ages of 1, 3, and 10 Myr. The two 
panels show the BCAH tracks for subsolar ([m/H] = -0.5; left) and solar (right) metallicities. 
In each of the panels for the BCAH and SDF tracks, the solid and dotted red lines indicate 
the interpolated tracks at our estimated dynamical masses of the primary and secondary of 
Haro l-Uc and their 1-a uncertainties, M x = 0.96±° f 8 M and M 2 = 0.33±g;g| M . The 
red circles show the location of the effective temperatures and absolute K-band magnitudes. 
The error bars in effective temperature represent ± 1 spectral subclass and those in absolute 
magnitude represent only the photometric uncertainties. The uncertainties in the distance 
modulus and are displayed by the vertical blue and purple error bars, respectively. 
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Fig. 9. — NIRSPEC spectrum in order 49 of Haro l-14c (solid line) compared with the 
spectral template HR 8085 (dotted line). The location of the Fe I lines are marked in the 
plot. 



